Iterative Landmark-Based Umbrella Sampling (ILBUS) Protocol for Sampling of Conformational Space of Biomolecules.

Computer simulations of biomolecules such as molecular dynamics often suffer from insufficient sampling. Due to limited computational resources, insufficient sampling prevents obtaining proper equilibrium distributions of observed properties. To deal with this problem, we proposed a simulation protocol for efficient resampling of collected off-equilibrium trajectories. These trajectories are utilized for the initial mapping of the conformational space, which is later properly resampled by the introduced Iterative Landmark-Based Umbrella Sampling (ILBUS) method. Reconstruction of static equilibrium properties is achieved by the multistate Bennett acceptance ratio (MBAR) method, which enables efficient use of simulated data. The ILBUS protocol is geometry-based and does not demand any additional collective variable or a dimensional-reduction technique. The only requirement is a set of suitably spaced reference conformations, which serve as landmarks in the mapped conformational space. Additionally, the ILBUS protocol encompasses an iterative process that optimizes the force constant used in the umbrella sampling simulation. Such tuning is an inherent feature of the protocol and does not need to be performed by the user in advance. Furthermore, even the simulations with suboptimal force constants can be used in estimates by MBAR. We demonstrate the feasibility and the performance of this approach in the study of the conformational landscape of the alanine dipeptide, met-enkephalin, and adenylate kinase.

Fusion of two unrelated protein domains in a chimera protein and its 3D prediction: Justification of the x-ray reference structures as a prediction benchmark.

Proteins are naturally formed by domains edging their functional and structural properties. A domain out of the context of an entire protein can retain its structure and to some extent also function on its own. These properties rationalize construction of artificial fusion multidomain proteins with unique combination of various functions. Information on the specific functional and structural characteristics of individual domains in the context of new artificial fusion proteins is inevitably encoded in sequential order of composing domains defining their mutual spatial positions. So the challenges in designing new proteins with new domain combinations lie dominantly in structure/function prediction and its context dependency. Despite the enormous body of publications on artificial fusion proteins, the task of their structure/function prediction is complex and nontrivial. The degree of spatial freedom facilitated by a linker between domains and their mutual orientation driven by noncovalent interactions is beyond a simple and straightforward methodology to predict their structure with reasonable accuracy. In the presented manuscript, we tested methodology using available modeling tools and computational methods. We show that the process and methodology of such prediction are not straightforward and must be done with care even when recently introduced AlphaFold II is used. We also addressed a question of benchmarking standards for prediction of multidomain protein structures-x-ray or Nuclear Magnetic Resonance experiments. On the study of six two-domain protein chimeras as well as their composing domains and their x-ray structures selected from PDB, we conclude that the major obstacle for justified prediction is inappropriate sampling of the conformational space by the explored methods. On the other hands, we can still address particular steps of the methodology and improve the process of chimera proteins prediction.

Modern and prebiotic amino acids support distinct structural profiles in proteins.

The earliest proteins had to rely on amino acids available on early Earth before the biosynthetic pathways for more complex amino acids evolved. In extant proteins, a significant fraction of the 'late' amino acids (such as Arg, Lys, His, Cys, Trp and Tyr) belong to essential catalytic and structure-stabilizing residues. How (or if) early proteins could sustain an early biosphere has been a major puzzle. Here, we analysed two combinatorial protein libraries representing proxies of the available sequence space at two different evolutionary stages. The first is composed of the entire alphabet of 20 amino acids while the second one consists of only 10 residues (ASDGLIPTEV) representing a consensus view of plausibly available amino acids through prebiotic chemistry. We show that compact conformations resistant to proteolysis are surprisingly similarly abundant in both libraries. In addition, the early alphabet proteins are inherently more soluble and refoldable, independent of the general Hsp70 chaperone activity. By contrast, chaperones significantly increase the otherwise poor solubility of the modern alphabet proteins suggesting their coevolution with the amino acid repertoire. Our work indicates that while both early and modern amino acids are predisposed to supporting protein structure, they do so with different biophysical properties and via different mechanisms.

TRPM5 Channel Binds Calcium-Binding Proteins Calmodulin and S100A1.

Melastatin transient receptor potential (TRPM) channels belong to one of the most significant subgroups of the transient receptor potential (TRP) channel family. Here, we studied the TRPM5 member, the receptor exposed to calcium-mediated activation, resulting in taste transduction. It is known that most TRP channels are highly modulated through interactions with extracellular and intracellular agents. The binding sites for these ligands are usually located at the intracellular N- and C-termini of the TRP channels, and they can demonstrate the character of an intrinsically disordered protein (IDP), which allows such a region to bind various types of molecules. We explored the N-termini of TRPM5 and found the intracellular regions for calcium-binding proteins (CBPs) the calmodulin (CaM) and calcium-binding protein S1 (S100A1) by in vitro binding assays. Furthermore, molecular docking and molecular dynamics simulations (MDs) of the discovered complexes confirmed their known common binding interface patterns and the uniqueness of the basic residues present in the TRPM binding regions for CaM/S100A1.

Amino Acid Interactions (INTAA) web server v2.0: a single service for computation of energetics and conservation in biomolecular 3D structures.

Interactions among amino acid residues are the principal contributor to the stability of the three-dimensional structure of a protein. The Amino Acid Interactions (INTAA) web server (https://bioinfo.uochb.cas.cz/INTAA/) has established itself as a unique computational resource, which enables users to calculate the contribution of individual residues in a biomolecular structure to its total energy using a molecular mechanical scoring function. In this update, we describe major additions to the web server which help solidify its position as a robust, comprehensive resource for biomolecular structure analysis. Importantly, a new continuum solvation model was introduced, allowing more accurate representation of electrostatic interactions in aqueous media. In addition, a low-overhead pipeline for the estimation of evolutionary conservation in protein chains has been added. New visualization options were introduced as well, allowing users to easily switch between and interrelate the energetic and evolutionary views of the investigated structures.

Intrinsically disordered protein domain of human ameloblastin in synthetic fusion with calmodulin increases calmodulin stability and modulates its function.

Constantly increasing attention to bioengineered proteins has led to the rapid development of new functional targets. Here we present the biophysical and functional characteristics of the newly designed CaM/AMBN-Ct fusion protein. The two-domain artificial target consists of calmodulin (CaM) and ameloblastin C-terminus (AMBN-Ct). CaM as a well-characterized calcium ions (Ca2+) binding protein offers plenty of options in terms of Ca2+ detection in biomedicine and biotechnologies. Highly negatively charged AMBN-Ct belongs to intrinsically disordered proteins (IDPs). CaM/AMBN-Ct was designed to open new ways of communication synergies between the domains with potential functional improvement. The character and function of CaM/AMBN-Ct were explored by biophysical and molecular modelling methods. Experimental studies have revealed increased stability and preserved CaM/AMBN-Ct function. The results of molecular dynamic simulations (MDs) outlined different interface patterns between the domains with potential allosteric communication within the fusion.

Characterization of AMBN I and II Isoforms and Study of Their Ca2+-Binding Properties.

Ameloblastin (Ambn) as an intrinsically disordered protein (IDP) stands for an important role in the formation of enamel-the hardest biomineralized tissue commonly formed in vertebrates. The human ameloblastin (AMBN) is expressed in two isoforms: full-length isoform I (AMBN ISO I) and isoform II (AMBN ISO II), which is about 15 amino acid residues shorter than AMBN ISO I. The significant feature of AMBN-its oligomerization ability-is enabled due to a specific sequence encoded by exon 5 present at the N-terminal part in both known isoforms. In this study, we characterized AMBN ISO I and AMBN ISO II by biochemical and biophysical methods to determine their common features and differences. We confirmed that both AMBN ISO I and AMBN ISO II form oligomers in in vitro conditions. Due to an important role of AMBN in biomineralization, we further addressed the calcium (Ca2+)-binding properties of AMBN ISO I and ISO II. The binding properties of AMBN to Ca2+ may explain the role of AMBN in biomineralization and more generally in Ca2+ homeostasis processes.

Toward Ab Initio Protein Folding: Inherent Secondary Structure Propensity of Short Peptides from the Bioinformatics and Quantum-Chemical Perspective.

By combining bioinformatics with quantum-chemical calculations, we attempt to address quantitatively some of the physical principles underlying protein folding. The former allowed us to identify tripeptide sequences in existing protein three-dimensional structures with a strong preference for either helical or extended structure. The selected representatives of pro-helical and pro-extended sequences were converted into "isolated" tripeptides-capped at N- and C-termini-and these were subjected to an extensive conformational sampling and geometry optimization (typically thousands to tens of thousands of conformers for each tripeptide). For each conformer, the QM(DFT-D3)/COSMO-RS free-energy value was then calculated, Gconf(solv). The Δ Gconf(solv) is expected to provide an objective, unbiased, and quantitatively accurate measure of the conformational preference of the particular tripeptide sequence. It has been shown that irrespective of the helical vs extended preferences of the selected tripeptide sequences in context of the protein, most of the low-energy conformers of isolated tripeptides prefer the R-helical structure. Nevertheless, pro-helical tripeptides show slightly stronger helix preference than their pro-extended counterparts. Furthermore, when the sampling is repeated in the presence of a partner tripeptide to mimic the situation in a ß-sheet, pro-extended tripeptides (exemplified by the VIV) show a larger free-energy benefit than pro-helical tripeptides (exemplified by the EAM). This effect is even more pronounced in a hydrophobic solvent, which mimics the less polar parts of a protein. This is in line with our bioinformatic results showing that the majority of pro-extended tripeptides are hydrophobic. The preference for a specific secondary structure by the studied tripeptides is thus governed by the plasticity to adopt to its environment. In addition, we show that most of the "naturally occurring" conformations of tripeptide sequences, i.e., those found in existing three-dimensional protein structures, are within ∼10 kcal·mol-1 from their global minima. In summary, our "ab initio" data suggest that complex protein structures may start to emerge already at the level of their small oligopeptidic units, which is in line with a hierarchical nature of protein folding.

AMBER and CHARMM Force Fields Inconsistently Portray the Microscopic Details of Phosphorylation.

Phosphorylation of serine, threonine, and tyrosine is one of the most frequently occurring and crucial post-translational modifications of proteins often associated with important structural and functional changes. We investigated the direct effect of phosphorylation on the intrinsic conformational preferences of amino acids as a potential trigger of larger structural events. We conducted a comparative study of force fields on terminally capped amino acids (dipeptides) as the simplest model for phosphorylation. Our bias-exchange metadynamics simulations revealed that all model dipeptides sampled a great heterogeneity of ensembles affected by introduction of mono- and dianionic phosphate groups. However, the detected changes in populations of backbone conformers and side-chain rotamers did not reveal a strong discriminatory shift in preferences, as could be anticipated for the bulky, charged phosphate group. Furthermore, the AMBER and CHARMM force fields provided inconsistent populations of individual conformers as well as net structural trends upon phosphorylation. Detailed analysis of ensembles revealed competition between hydration and formation of internal hydrogen bonds involving amide hydrogens and the phosphate group. The observed difference in hydration free energy and potential for hydrogen bonding in individual force fields could be attributed to the different partial atomic charges used in each force field and, hence, the different parametrization strategies. Nevertheless, conformational propensities and net structural changes upon phosphorylation are difficult to extract from experimental measurements, and existing experimental data provide limited guidance for force field assessment and further development.

Sequence Versus Composition: What Prescribes IDP Biophysical Properties?

Intrinsically disordered proteins (IDPs) represent a distinct class of proteins and are distinguished from globular proteins by conformational plasticity, high evolvability and a broad functional repertoire. Some of their properties are reminiscent of early proteins, but their abundance in eukaryotes, functional properties and compositional bias suggest that IDPs appeared at later evolutionary stages. The spectrum of IDP properties and their determinants are still not well defined. This study compares rudimentary physicochemical properties of IDPs and globular proteins using bioinformatic analysis on the level of their native sequences and random sequence permutations, addressing the contributions of composition versus sequence as determinants of the properties. IDPs have, on average, lower predicted secondary structure contents and aggregation propensities and biased amino acid compositions. However, our study shows that IDPs exhibit a broad range of these properties. Induced fold IDPs exhibit very similar compositions and secondary structure/aggregation propensities to globular proteins, and can be distinguished from unfoldable IDPs based on analysis of these sequence properties. While amino acid composition seems to be a major determinant of aggregation and secondary structure propensities, sequence randomization does not result in dramatic changes to these properties, but for both IDPs and globular proteins seems to fine-tune the tradeoff between folding and aggregation.

Widespread evolutionary crosstalk among protein domains in the context of multi-domain proteins.

Domains are distinct units within proteins that typically can fold independently into recognizable three-dimensional structures to facilitate their functions. The structural and functional independence of protein domains is reflected by their apparent modularity in the context of multi-domain proteins. In this work, we examined the coupling of evolution of domain sequences co-occurring within multi-domain proteins to see if it proceeds independently, or in a coordinated manner. We used continuous information theory measures to assess the extent of correlated mutations among domains in multi-domain proteins from organisms across the tree of life. In all multi-domain architectures we examined, domains co-occurring within protein sequences had to some degree undergone concerted evolution. This finding challenges the notion of complete modularity and independence of protein domains, providing new perspective on the evolution of protein sequence and function.

[Life threatening manifestations of lupus and antiphospholipid syndrome in internal medicine]. / Zivot ohrozující projevy systémového lupusu a antifosfolipidového syndromu ve vnitrním lékarství.

The clinical picture of systemic lupus and antiphospholipid syndrome is remarkably varied and disease manifestations are commonly very heterogeneous. Relatively often both diseases are associated with severe, acute and life threatening manifestations, which places demands on the knowledge of differential diagnostics and experience of the physicians. This article deals with the serious and mostly acute impairment of cardiovascular, respiratory, renal, gastrointestinal, hematopoietic or nervous systems, briefly discusses the acute pregnancy complication and summarizes the basic therapeutic option. It emphasizes the role of both, sometimes inseparable, diseases in differential diagnosis of acute symptoms in internal medicine.Key words: clinical symptoms - diagnostics - live threatening manifestations - lupus erythematosus - systemic antiphospholipid syndrome - therapy.

Cardiac amyloidosis: from clinical suspicion to morphological diagnosis.

Amyloidosis is a heterogeneous group of diseases characterised by extracellular accumulation of amyloid in various tissues and organs of the body, leading to alteration and destruction of tissues. Heart involvement is the most important prognostic factor in patients with systemic amyloidosis and the diagnosis and typing of amyloid must be made properly. The clinical picture shows congestive heart failure with predominant right-sided heart failure symptoms in fully developed disease, various types of arrhythmias and characteristic electrocardiography and echocardiography findings. Blood and urine monoclonal protein studies and cardiac biomarkers belong to the spectrum of standard laboratory examinations. Cardiac cardiomyopathy is connected with amyloid based on immunoglobulin light chains, serum amyloid A, transthyretin, atrial natriuretic factor or apolipoprotein A1. In the routine diagnostic algorithm, biopsy specimens are examined using special histological staining, immunohistochemistry and immunofluorescence; proteomic analysis is only performed in specialised centres.

Random protein sequences can form defined secondary structures and are well-tolerated in vivo.

The protein sequences found in nature represent a tiny fraction of the potential sequences that could be constructed from the 20-amino-acid alphabet. To help define the properties that shaped proteins to stand out from the space of possible alternatives, we conducted a systematic computational and experimental exploration of random (unevolved) sequences in comparison with biological proteins. In our study, combinations of secondary structure, disorder, and aggregation predictions are accompanied by experimental characterization of selected proteins. We found that the overall secondary structure and physicochemical properties of random and biological sequences are very similar. Moreover, random sequences can be well-tolerated by living cells. Contrary to early hypotheses about the toxicity of random and disordered proteins, we found that random sequences with high disorder have low aggregation propensity (unlike random sequences with high structural content) and were particularly well-tolerated. This direct structure content/aggregation propensity dependence differentiates random and biological proteins. Our study indicates that while random sequences can be both structured and disordered, the properties of the latter make them better suited as progenitors (in both in vivo and in vitro settings) for further evolution of complex, soluble, three-dimensional scaffolds that can perform specific biochemical tasks.

Focused conformational sampling in proteins.

A detailed understanding of the conformational dynamics of biological molecules is difficult to obtain by experimental techniques due to resolution limitations in both time and space. Computer simulations avoid these in theory but are often too short to sample rare events reliably. Here we show that the progress index-guided sampling (PIGS) protocol can be used to enhance the sampling of rare events in selected parts of biomolecules without perturbing the remainder of the system. The method is very easy to use as it only requires as essential input a set of several features representing the parts of interest sufficiently. In this feature space, new states are discovered by spontaneous fluctuations alone and in unsupervised fashion. Because there are no energetic biases acting on phase space variables or projections thereof, the trajectories PIGS generates can be analyzed directly in the framework of transition networks. We demonstrate the possibility and usefulness of such focused explorations of biomolecules with two loops that are part of the binding sites of bromodomains, a family of epigenetic "reader" modules. This real-life application uncovers states that are structurally and kinetically far away from the initial crystallographic structures and are also metastable. Representative conformations are intended to be used in future high-throughput virtual screening campaigns.

Amyloid ß Fibril Elongation by Monomers Involves Disorder at the Tip.

The growth of amyloid fibrils from Aß1-42 peptide, one of the key pathogenic players in Alzheimer's disease, is believed to follow a nucleation-elongation mechanism. Fibril elongation is often described as a "dock-lock" procedure, where a disordered monomer adsorbs to an existing fibril in a relatively fast process (docking), followed by a slower conformational transition toward the ordered state of the template (locking). Here, we use molecular dynamics simulations of an ordered pentamer of Aß42 at fully atomistic resolution, which includes solvent, to characterize the elongation process. We construct a Markov state model from an ensemble of short trajectories generated by an advanced sampling algorithm that efficiently diversifies a subset of the system without any bias forces. This subset corresponds to selected dihedral angles of the peptide chain at the fibril tip favored to be the fast growing one experimentally. From the network model, we extract distinct locking pathways covering time scales in the high microsecond regime. Slow steps are associated with the exchange of hydrophobic contacts, between nonnative and native intermolecular contacts as well as between intra- and intermolecular ones. The N-terminal segments, which are disordered in fibrils and typically considered inert, are able to shield the lateral interfaces of the pentamer. We conclude by discussing our findings in the context of a refined dock-lock model of Aß fibril elongation, which involves structural disorder for more than one monomer at the growing tip.

Amino Acid Interaction (INTAA) web server.

Large biomolecules-proteins and nucleic acids-are composed of building blocks which define their identity, properties and binding capabilities. In order to shed light on the energetic side of interactions of amino acids between themselves and with deoxyribonucleotides, we present the Amino Acid Interaction web server (http://bioinfo.uochb.cas.cz/INTAA/). INTAA offers the calculation of the residue Interaction Energy Matrix for any protein structure (deposited in Protein Data Bank or submitted by the user) and a comprehensive analysis of the interfaces in protein-DNA complexes. The Interaction Energy Matrix web application aims to identify key residues within protein structures which contribute significantly to the stability of the protein. The application provides an interactive user interface enhanced by 3D structure viewer for efficient visualization of pairwise and net interaction energies of individual amino acids, side chains and backbones. The protein-DNA interaction analysis part of the web server allows the user to view the relative abundance of various configurations of amino acid-deoxyribonucleotide pairs found at the protein-DNA interface and the interaction energies corresponding to these configurations calculated using a molecular mechanical force field. The effects of the sugar-phosphate moiety and of the dielectric properties of the solvent on the interaction energies can be studied for the various configurations.

Emergency situations in rheumatology with a focus on systemic autoimmune diseases.

BACKGROUND AND AIM: Rheumatic diseases are commonly considered chronic conditions. However, acute manifestations can be very severe and represent a diagnostic problem. Examples are systemic lupus erythematosus with acute flare, glomerulonephritis, CNS disorders and catastrophic antiphospholipid syndrome, scleroderma with interstitial lung disease, pulmonary hypertension and renal crisis and polyangiitis with alveolar haemorhage and acute respiratory failure. This aim of this paper is to overview emergency situations which can be encountered in the care of patients with autoimmune systemic diseases and vasculitides. METHODS: A Pubmed search for both original and review articles, recent textbooks and current guidelines related to rheumatic diseases with possible acute situations were included in this review article. Relevant image documentation was obtained at the site over the past several years of observation. CONCLUSIONS: This paper provides an overview of facts and emergency situations which can be encountered in the care of patients with autoimmune systemic diseases and vasculitides. It is directed at clinicians working in intensive care. It provides a differential diagnostic overview and information which is rare and commonly underestimated.

Effect of TFE on the Helical Content of AK17 and HAL-1 Peptides: Theoretical Insights into the Mechanism of Helix Stabilization.

Fluorinated alcohols such as 2,2,2-trifluoroethanol (TFE) are among the most frequently used cosolvents in experiment studies of peptides. They have significant effects on secondary structure and a particularly strong promotion of α-helix is induced by TFE. In this study we validated recently proposed force field parameters for TFE in molecular dynamics simulations with two model peptides-alanine-rich AK-17 and antimicrobial peptide halictine-1 (HAL-1). In the case of HAL-1, we characterized the effect of TFE on this peptide experimentally by ECD spectroscopy. Our TFE model in question reproduced the helix-promoting effect of TFE and provided insight into the mechanisms of TFE action on peptides. Our simulations confirmed the preferential interaction of TFE molecules with α-helices, although the TFE molecules accumulate in the vicinity of the peptides in various conformations. Moreover, we observed a significant effect of TFE on the thermodynamics of the helix-coil transition and a change in local conformational preferences in the unfolded (coil) state induced by TFE. In addition, our simulation-based analysis suggests that different mechanisms participate in helix stabilization in both model peptides in water and TFE solution. Our results thus support the picture of complex TFE action on peptides that is further diversified by the identity and intrinsic properties of the peptide.

Insights into Unfolded Proteins from the Intrinsic ϕ/ψ Propensities of the AAXAA Host-Guest Series.

Various host-guest peptide series are used by experimentalists as reference conformational states. One such use is as a baseline for random-coil NMR chemical shifts. Comparison to this random-coil baseline, through secondary chemical shifts, is used to infer protein secondary structure. The use of these random-coil data sets rests on the perception that the reference chemical shifts arise from states where there is little or no conformational bias. However, there is growing evidence that the conformational composition of natively and nonnatively unfolded proteins fail to approach anything that can be construed as random coil. Here, we use molecular dynamics simulations of an alanine-based host-guest peptide series (AAXAA) as a model of unfolded and denatured states to examine the intrinsic propensities of the amino acids. We produced ensembles that are in good agreement with the experimental NMR chemical shifts and confirm that the sampling of the 20 natural amino acids in this peptide series is be far from random. Preferences toward certain regions of conformational space were both present and dependent upon the environment when compared under conditions typically used to denature proteins, i.e., thermal and chemical denaturation. Moreover, the simulations allowed us to examine the conformational makeup of the underlying ensembles giving rise to the ensemble-averaged chemical shifts. We present these data as an intrinsic backbone propensity library that forms part of our Structural Library of Intrinsic Residue Propensities to inform model building, to aid in interpretation of experiment, and for structure prediction of natively and nonnatively unfolded states.